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Abstract 

The statistics of persistent events, recently introduced in the context of phase 
ordering dynamics, is investigated in the case of the one-dimensional lattice ran- 
dom walk in discrete time. We determine the survival probability of the random 
walker in the presence of an obstacle moving ballistically with velocity v, i.e., the 
probability that the random walker remains up to time n on the left of the obstacle. 
Three regimes are to be considered for the long-time behavior of this probability, 
according to the sign of the difference between v and the drift velocity V of the 
random walker. In one of these regimes (v > V), the survival probability has a 
non-trivial limit at long times, which is discontinuous at all rational values of v. 
An algebraic approach allows us to compute these discontinuities, as well as sev- 
eral related quantities. The mathematical structure underlying the solvability of 
this model combines elementary number theory, algebraic functions, and algebraic 
curves defined over the rationals. 
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1 Introduction 



Consider an asymmetric binomial random walk on a one-dimensional lattice. In units 
of the lattice spacing, the steps e m performed by the walker at integer times t = m are 
independent identically distributed random variables, with the binary law 

{—1 with probability p, 
(1.1) 
+1 with probability q = 1 — p. 

This paper is devoted to the analysis of the probability F(n,v) that the walker re- 
mains, up to time n, on the left of an obstacle moving ballistically with velocity v. In a 
pictorial language, F(n, v) is the probability that a sober man walking at constant speed 
v, and a drunken sailor stepping forward and backward erratically, both leaving a pub 
at some initial time, do not meet up to time n. Alternatively, F(n,v) can be viewed as 
the survival probability of the walker in the presence of the obstacle. Surprisingly, as 
we shall see, this probability is highly non-trivial, especially as far as its w-dependence is 
concerned. 

The quantity F(n, v) also represents the probability that the path of the random 
walker in the (t, x)-plane remains, up to the integer time t = n, on the left of the straight 
wall x = vt. Figures 1 to 3 provide illustrations of these definitions. 




Figure 1: Configuration space of the problem. 

This investigation was motivated by recent work on the statistics of persistent events 
in nonequilibrium statistical-mechanical systems undergoing phase ordering [lj. Persis- 
tent events are defined by a constraint on the past history of the system, to be satisfied 
up to time t, i.e., for all previous times < t' < t. This concept can be simply illustrated 
on the example of a chain of Ising spins, starting from a random initial condition, and 
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evolving under Glauber or heat-bath dynamics at zero temperature O. The persistence 



probability R(t) is defined as the fraction of space which remained in the same phase 
up to time t, i.e., in the present the fraction of spins which did not flip up to 

time t 0. For long times the persistence probability decays as 



where 6 is the persistence exponent. 




Figure 2: A walk remaining on the left of the wall, thus contributing to the survival 
probability. 




Figure 3: A walk crossing the wall, thus not contributing to the survival probability. 

The persistence exponent is actually only one member of a continuous family of 
exponents defined as follows Denoting by a(t) the spin at a given site, consider the 
local mean magnetization 



R(t) ~ t 



e 



(1.2) 




(1.3) 
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This quantity is simply related to the fraction of time spent by the spin in the positive 
direction. Consider now the probability of persistent large deviations R(t, x) above the 
level x (with — 1 < x < 1), defined as the probability that M(t) remained greater than x, 
for all times < t' < t. For the Glauber-Ising chain at zero temperature, this quantity 
is observed to decay algebraically at large times, with an exponent 8(x) continuously 
varying with x When x — 1, the usual persistence probability is recovered, so that 
#(1) = 6. The exponent 9{x) appears as a first-passage exponent associated to persistent 
large deviations of M(t) U g, g. 

The interpretation of the stochastic process a(t) as the steps of a fictitious random 
walker naturally led us to the present study, where, for simplicity, we consider the con- 
cepts introduced above in the simplest case of a binomial random walk. In the present 
work, the analogue of M(t) is the mean velocity V n of the walker at time n, 

V n = ^, (1.4) 

n 

where 

n 

%n = &m (1-5) 

m=l 

is the position of the walker at time n, while F(n, v) plays the role of the probability of 
persistent large deviations R(t,x). Note a slight difference in the above definitions. In 
this work, F(n,v) is the probability that the mean speed was always less than v, while 
in refs. Jj], || R(t, x) is the probability that the mean magnetization was always above 
x. This difference is harmless, as it just amounts to interchanging the probabilities p and 
q and changing the slope v into its opposite. 

More precise definitions of the quantities considered in this work are as follows. First, 
the 'one-time' distribution functions P ± (n, v) of the mean velocity V n at time n are 
defined for a given slope v as 

P~(n,v) = Prob{V n < v} = Prob{x n < nv}, 

P + (n,v) = Prob{\4 < v} — Prob{x n < nv}, (1.6) 

Then the survival probabilities F ± (n, v) up to time n are defined as 

F~(n, v) = ProbjVm < v for m = 1, . . . ,n} = Prob{x m < mv for m = 1, . . . , n}, 
F + (n, v) = Prob{V^i < v for m = 1, . . . , n} = Prob{x m < mv for m = 1, . . . , n}. 

(1.7) 

Two definitions are needed for each quantity, since V n is a discrete random variable, so 
that P ± (n, v) and F ± (n, v) are in general discontinuous at any rational value of the slope 
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v. Throughout this paper, quantities with the + and the — superscript are distinct from 
each other only when the slope v is rational. Equations involving P ± (n,v), or P(n,v) 
for short, are meant to hold separately for P + (n,v) and P~(n,v). 

We find three regimes for the long-time behavior of P ± (n, v) and F ± (n, v), according 
to the sign of the relative velocity of the obstacle with respect to the drift velocity of the 
random walker 

V = lim V n = (e) = 1 - 2p, (1.8) 

n^oo 

with the following results. 

• Large-deviation regime (v < V): 

Both P ± (n,v) and F (n, v) decay exponentially in time, with a common entropy 
function S(v), but different power- law prefactors, namely 

P ± (n, v) « a±(v) n~ 1/2 e~ nS{v \ (1.9) 
F ± (n,v) ^b±(v)n- 3/2 e- nS{v \ (1.10) 

• Marginal regime (v — V): 
We have 

P ± (n,v) -> ^ (1.11) 

(1-12) 
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The limit value 1/2 of P ± (n,v) is a consequence of the central limit theorem. 
The inverse-square-root decay of F ± (n, v) is a known property of one-dimensional 
random walks || [l], HI • 

• Convergent regime (v > V): 

P ± (n, v) converge exponentially to one, while P ± (n, v) admit non-trivial limits: 

1 - P ± (n, v) w a± (w)n" 1/2 e - nS ^\ (1.13) 
^(n,!;) -> ^(v). (1.14) 

These results were partially announced in ref. M. The asymptotic behavior of the 
one-time distribution functions P ± (n,v) are well-known properties of random walks ||, 
while the long-time behavior of the survival probabilities F ± (n, v) in these three regimes 
is the subject of the present study. To obtain this long-time behavior, we shall proceed 
in two steps. 
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A first level of description is provided in sec. 0. We use a combinatorial result, 
originally due to Sparre Andersen, and proved in Appendix A, in order to derive the 
form of the estimates ( |1.10| ), ( |1.12| ), and (fTT4j). The prefactors b^{v) and C (v), as well 
as the limit survival probabilities F ± (v), are respectively given by eqs. ( |2.22| ), ( |2.28| ), 
and ( |2.33|) . Such expressions hold for random walks with any distribution of steps e m , 
whose first two moments are finite. However they are only formal, in the sense that they 
do not lead to closed-form results in general. 

A second, more refined level of description is provided by the rest of the paper, de- 
voted to a detailed study of the quantities b^(v), C ± (f), and F ± (v). This investigation 
will allow us to unravel an unexpected richness in the statistics of persistent events in 
the present situation of the binomial random walk on the lattice of integers. More pre- 
cisely, sec. |3| contains basic concepts and results, emphasizing the central importance of 
sequences of positive integers Af. Sec. |]is devoted to the exposition of several methods 
(continuity of the path, probability flow, duality symmetry), which allow the determina- 
tion of the integers and related quantities, in the case of an arbitrary slope v. The 
case of a rational slope v is dealt with in sec. [5], where algebraic functions play a central 
role. Sec. ^| contains an investigation of the critical behavior of the survival probabilities, 
obtaining thus predictions in the three regimes. Further results on several specific situ- 
ations, including the slopes of the form v = 1 — 2/N or v = — 1 + 2/N, are derived in 
sec. |7|, while sec. || contains a discussion. 



2 Sparre Andersen formalism 

As announced in the Introduction, this section is devoted to a first level of description of 
the one-time distribution functions P ± (n, v) and of the survival probabilities F ± (n, v). 

We shall make an extensive use of a remarkable combinatorial result, originally due 
to Sparre Andersen 0. The possibility of recasting the computation of F ± (n, v) as a 
Sparre Andersen problem was uncovered in ref. Q. The result of Sparre Andersen is 
best expressed as an identity between the generating series of the P ± (n,v) and of the 
F ± (n, v ), namely 

= £ F^n^z- = exp ^hA z A . (2.1) 

n>0 \n>l 71 j 

This result can be found e.g. in the book by Feller ||. Since eq. (|2.1| ) lacks a clear 
intuitive meaning, we give an elementary and self-contained combinatorial proof of it in 
Appendix A. 
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Taking the logarithmic derivative of both sides of eq. ( j2.1|) , we obtain linear recursion 
relations for F ± (n, v): 

n-l 

nF ± (n, v ) = p± ( n - m > v)F ± (m, v), (2.2) 

m=0 

with F ± (0,f) = 1, hence 



F*(2,v) = i((P ± (l,^)) 2 + P ± (2 ! . 



F ± (3,v) = ^((F^M)) 8 + 3P ± (l,v)P ± (2,v) + 2P ± (3,v)) , (2.3) 

and so on. 

In the present situation, the one-time distribution functions P ± (ra, v) can be evaluated 
in closed form as follows. We introduce the notations 

1 — v 1 + v 

Pc = — 2~> q c = l-p c = — • (2-4) 

Let k (respectively, k) be the number of steps to the left (respectively, to the right), 
i.e., the number of e m equal to —1 (respectively, to +1), among the first n steps. The 
probability distribution of the integer k — 0, . . . , n reads 

Pn, k = (fyp k q n - k . (2-5) 
The position x n of the particle at time n is related to k and k by 

n = k + k, x n = k — k = n — 2k, k = — 7T~~i k = — ~L ~ • (2-6) 

The condition V n < v (respectively, V n < v) is equivalent to k > np c (respectively, 
k > np c ), hence the result 

n 

P ± (n,v)= Pn,k, (2.7) 

fc=Int ± (np c )+l 

where we have defined the following two integer-part and fractional-part functions: 

, , , f < Frac"(x) < 1, 

x = Int (x) + Frac (x), with Int (x) integer and < , (2-8) 

( < Frac (x) < 1. 

The random walk possesses a duality symmetry, defined by interchanging the prob- 
abilities p and q, and simultaneously changing the slope v into its opposite —v, i.e., 
interchanging p c and q c . This amounts to considering the survival probability in the 
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complementary domain, i.e., on the right of the wall shown in Figure 1. Let us make 
the dependence of quantities on p explicit for a while, indicating the value of p after 
a semi-colon. Our conventions imply that the one-time distribution functions pertain- 
ing to the left of the wall and to the right of the wall read respectively P ± (n,v;p) and 
1 — P T (n, —v; q). These quantities are clearly equal: 

P ± (n,v;p) + P T {n, -v;q) = 1. (2.9) 

As a consequence, the generating series f ± (z 1 v) obeys the following identity: 

f±( z , v; p)f*( z ,- v;q ) = J—. (2.10) 



More spectacular consequences of this duality symmetry will be investigated in sees. JO 
and |5.3| , by means of more powerful techniques. 

We now turn to the investigation of the asymptotic behavior of the survival proba- 
bilities. As mentioned above, three regimes are to be considered, according to the sign 
of the relative velocity of the obstacle with respect to the drift velocity of the random 
walker [see eq. ( |1.8| )1, namely 

v-V = v-{l-2p) = 2(p-p c ). (2.11) 
2.1 Large-deviation regime: v < V or p < p c 

This situation is typical of large deviations. The chance for the mean velocity of the 
random walker to deviate from its average is exponentially decreasing with time. More 
precisely, a careful treatment of the sums in eq. ( |2.7| ), using eq. ( |2.5|) and Stirling's 
formula, leads to 

P ± {n, v) « a±(v) n- 1 ' 2 e~ nS{v) (v < V) (2.12) 



see eq. (|1.9f) ]. The associated entropy function (or large-deviation function) reads 



S (v) = p c In — + q c In — 
p q 

= H (i ^ )ini ^ +(i+ '" )ini ^)- (2j3) 

This function is positive, and it vanishes quadratically as p — > p c , i.e., V — > v, in 
agreement with the central limit theorem, according to 
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In eq. (|2.12|) , the prefactors a^(v) = a ± (np c ) are v -dependent periodic functions of np c 
with unit period, which can be determined explicitly: 



1-Frac±(x) 



a ± {x) = \P^J (2.15) 

(2™Wi-M 

V Pel J 

with the definition (|2.8|) . The amplitude functions a ± (x) can be alternatively expanded 
as Fourier series: 

oo 



a (x 



£ afe 2wUx , (2.16) 



e=-oo 



so that eq. (|2.12| ) can be recast as 

oo 

P ± {n,v) ^n- 1 ' 2 afe( 2niEpc - s{v) ) n . (2.17) 
e=-oo 

We remind for further reference the following general result, known as a Tauberian 
theorem ||. If positive numbers, with asymptotic behavior 

c n ^an^(z c y n (n>l), (2.18) 

then the power series 

f(z) = J2 Cnz n (2-19) 

n>0 

defines an analytic function whose nearest singularity is at z = z c , where it has a power- 
law singular part of the form 

/ sg ( Z )«ar( 7 + l) (l--) (z-* c -0), (2.20) 
\ z c / 

T denoting Euler's gamma function. The reciprocal property, namely that eq. ( |2.20| ) 
implies eq. (|2.18|) , holds if the sequence c n is assumed to be smooth enough. 

As a consequence of eq. ( p.20|) , the estimate ( p.!7|) implies that the series f ± (z,v) of 
eq. Q2.1| ) have square-root branch points at z^ = exp(S(v) — 2ni£p c ), where £ runs over 
the integers. Using eq. (|2.18|) , we obtain the estimate 

F ± (n,v) « b±{v)n-* /2 e- nS{v) (v < V), (2.21) 

with b^(v) = b ± (np c ), where b ± (x) are again periodic functions, with unit period, whose 
Fourier coefficients are related to those of a ± (a;) by 



bf = aff ± (z e ). (2.22) 
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The survival probabilities F ± (n,v) thus decay exponentially in the large-deviation 
regime, with the same entropy function as the one-time distribution functions P ± (n,v), 
but with a different power of n [see eqs. QL9D , ( |1.10| )1, multiplied by non-trivial periodic 
prefactors b ± (np c ). 
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Figure 4: Periodic amplitudes a(x) and b(x), for p = 0.3 and v equal to the golden 
slope. 

Figure 4 shows a plot of the periodic amplitudes a(x) and b(x), for p = 0.3, and for 
the 'golden slope' v, such that p c is the inverse golden mean, the most typical irrational: 

V$- 1 



Pc 



0.618034, 



2-V5 



-0.236068. 



(2.23) 



The data for P(n, v) have been obtained by means of eq. (|2.7| ), and for F(n, v) by solving 
the recursion relation fl2.2p . The amplitude a(x) is found to be a pure exponential, in 
quantitative agreement with the prediction ( |2.15|) . The amplitude b(x) exhibits a rich 
structure, with discontinuities at all positive integer multiples of p c modulo 1. The 
occurrence of these discontinuities can be explained as follows. Eq. ( |2.2j ) implies that 

nF ± (n, v) = P ± (n, v) + F ± (l, v)P ± {n -l,v) + F ± (2, v)P ± {n - 2, v) + - • • (2.24) 

is a linear combination of P ± (n,v), P ± (n — l,v), P ± (n — 2,v), and so on, with n- 
independent coefficients. If we replace in eq. Q2.24Q P ± (n, v), P ± {n— 1, v), P ± (n—2, v), by 
their asymptotic form ( |2.12 ), with a„ (u) = a ± (np c ), we formally recover the result ( |2.2 
with 



b ± (x) =a ± (x) + F ± (l, v) e s{v) a ± (x - p c ) + F ± (2, v) e 2S{v) a ± (x - 2p c ) + • 



(2.25) 
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Since the amplitude functions a ± (x), given by eq. (|2.15|) , are discontinuous at x = 



modulo 1, eq. (|2.25|) shows that the amplitude functions b ± (x) are discontinuous at 
all positive integer multiples of p c modulo 1, confirming thus the observation made on 
Figure 4. 

2.2 Marginal regime: v = V or p = p c 

In this regime, the drift velocity V of the random walker is equal to the velocity v of the 
obstacle. The one-time distribution functions P ± (n,v) are expected to go to 1/2, as a 
consequence of the central limit theorem [see eq. (|1.11|) 1. 

A careful treatment of eq. (|2.7|) leads to the more accurate estimate 



P ± (n, v) w - + (^Frac ± (np c ) + Pc ^ (2nnp c q c y 1/2 (n > 1). (2.26) 



The Sparre Andersen relation ( |2.1| ) then leads to the behavior 



with 



C>H ex P (g f± '"-';» - . (2.28) 



Eq. fl2.18| ) in turn leads to 



F±fav)*0L (v = V). (2-29) 

The survival probabilities F ± (n, v) thus decay according to a universal inverse-square- 
root law in the marginal already mentioned in eq. ( |1.12| ). The associated ampli- 



tudes C ± (v ) are non-trivial functions of v, or equivalently of p c , to which we shall come 
back in sec. ^ The duality symmetry ( 2.10|) implies the relation 



C ± (w)C T (-t;) = 1. (2.30) 
2.3 Convergent regime: v > V or p > p c 

In this regime, the drift velocity V of the random walker belongs to the domain in 
which the mean velocities V n are constrained. These constraints are therefore less and 
less stringent as time goes on. Indeed P ± (n,v) go to unity exponentially fast, with 
1 — P ± (n,v) given by a large- deviation expression [see eq. ( 1.13| )]. 
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The Sparre Andersen relation ( |2.1|) now leads to the estimate 

F ± (v) 

^yrj (2.31) 

In other words, the survival probabilities converge toward non-trivial limits for infinite 
times: 

F ± (v) = lim F ± (n,v), (2.32) 

n— >oo 

hence the word 'convergent regime'. The limit survival probabilities have the formal 
expressions 

F ± (v) = exp (- Xj 1 ~ j (v>V). (2.33) 

The limit probabilities F (v ) have in general a very rich dependence on v and p, which 
is not apparent in the formal expressions (|2.33|) , and which will be investigated in the 



following sections. For the time being, we want to underline the existence of a universal 
relationship between the p-dependence of F ± (v) as p — > p c + and the amplitudes C ± (f ) 
corresponding to the marginal regime (p = p c ). For p — p c small and positive, the one- 
time distribution functions P ± (n,v) can be approximated by their values at p = p c as 
long as n(p — p c ) <C 1, while they can be estimated from the central limit theorem for 
n{p-p c ) ~ 1: 

' n(p-p c y 
, 2 3/2 p c q c 

where erf is the error function. By inserting the above estimate into eq. fl2.33|) , and 



P ± (n,v) ps - 



1 + erf 



(2.34) 



splitting the sum over n into three sums, respectively corresponding to the ranges 1 < 
n < n , no < n < ni, and n > n\, with 1 < % C < n>i, we obtain after some 

manipulations 



/ 2 \ 1/2 
*{v) « C^) (p-p c ) 



The survival probabilities F ± (f ) thus vanish linearly as p — > p c , with amplitudes propor- 
tional to the constants C ± (v), characteristic of the fall off ( |2.29| ) in the marginal regime 

(P=Pc)- 

The quantities introduced to describe the asymptotic long-time behavior of the sur- 
vival probabilities F ± (n, v) in the three different regimes, namely the periodic functions 
6 ± (x, v), the constants C ± (v), and the limit survival probabilities F ± (v), bear a highly 
non-trivial dependence on the parameters p and v, which is hidden in their formal ex- 
pressions (|2.22|) , ( 2.28 ), and ( [2.33|) . The rest of this paper is devoted to an investigation 
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of all these quantities, including their dependence on p and v. Throughout the following, 
when the context permits, we shall suppress the explicit dependence on v, writing e.g. 
F ± for F ± (v). 



3 Basic concepts 

3.1 The survival probability as a distribution function 

Define the maximal velocity V m3X of a given (infinite) random walk as the supremum of 
all the instantaneous mean velocities V n : 

V max = swpV n . (3.1) 

n 

This quantity is a random variable, as it depends on the walk under consideration. 
The survival probabilities F ± (w) can be interpreted as the distribution functions of this 
random variable, for a fixed value of p. We have indeed 

F~{v) = Prob{V max < v}, F+(v) = Prob{Knax < v}. (3.2) 

With probability one the random variable Vm ax is rational, and it lies in the range 
V < Vmax < 1- These properties originate in the following two facts. The V n converge 
with probability one to the drift velocity V, and the differences V n — V take both signs. 
As a consequence, with probability one there is some finite n such that V max = V n . We 
have therefore 

F ± {y) =0 (v < V), F ± (v) = 1 (v > 1). (3.3) 



The lower edge (v = V) corresponds to the marginal case studied in sec. |2.2| , while the 
upper edge {y = 1) will be investigated in sec. [7.1[ 

Whenever the slope v is irrational, we have F + (v) = F~(v), and the distribution 
function is continuous. On the contrary, when v is rational, quantities with superscripts + 
and — are different from each other in general. The discontinuity n(u) of the distribution 
function is nothing but the probability that the maximal velocity V max assumes the value 
v: 

U(v) = F + (v) - F-(v) = Prob{V max = v}, (3.4) 

and we have 

F + (v) = lim F(w), F~(v) = lim F(w). (3.5) 
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The above quantities can be related as follows. Let Q(v) be the probability that the 
random walker makes its first step to the left and then touches the wall at least once: 

Q{v) = Prob{V m < v for m = 1, . . . , n — 1 and V n = v}. (3.6) 

n>2 

This quantity is non-zero for any rational v in the range — 1 < v < 1. We can split any 
walk contributing to Q(v) into two independent walks: a finite walk from the origin to 
the first coincidence time n such that V n = v , and the rest, which is an infinite walk. 
Consider a walk contributing to F + (v). It either contributes to F~(v) or to 11(f). In the 
second case, we have V max = v, and with probability one there is some n such that V n — v 
for the first time. The complete walk splits into two independent walks: a finite walk from 
the origin to the first coincidence, and the rest, which is an infinite walk contributing to 
F + (v). So, in order to go from Q{y) to the difference F + {v) — F~(v), one just replaces 
the arbitrary infinite part by a walk contributing to F + (v ). In other words, we have 

U(v) = F + (v) - F~(v) = Q(v)F + (v), (3.7) 

with this quantity vanishing for v < V and being positive for v > V, or else 

F~(v) = (l-Q(v))F + (v). (3.8) 

The case of a rational slope v will be investigated in sec. £|. It turns out that some 
general combinatorial statements can be conveniently presented by keeping v arbitrary. 
This is the purpose of the rest of sec. |3|. 



3.2 Truncated Pascal triangle 

The oriented lattice edges crossing the wall from left to right will play an important role 
in what follows. Indeed the walks which contribute to the survival probability, i.e., those 
which remain on the left forever, are those which start to the left and never pass through 
a crossing edge. 

We assume for a while that v is irrational. As shown in Figure 5, we index the crossing 



edges by the integer k > 1, i.e., by the number of steps to the left, defined in eq. (|2.6| ). 
The endpoint of the k-th crossing edge has coordinates (n,x) = (nk,rik — 2k), where 
is defined by the inequalities k/p c < nu < k/p c + 1, or equivalently 

n k = l + lnt{k/p c ). (3.9) 

When v is rational, we have to be more careful and define two sequences. The edges 
corresponding to F~(v) may end but not start on the wall, while those corresponding to 
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F + (v) may start but not end on the wall, hence the prescriptions 



l + Int + (fc/p c ), 




l + Int-(A;/p c ). 



(3.10) 



Hence, if v is rational, any construction involving the sequence of crossing edges will 
give two different outputs in general, depending on which family of crossing edges is 
considered. In this case, we may, and shall sometimes, start the sequence of crossing 
edges with k = 0. 



Figure 5: Labeling of crossing edges by an integer k = 1, 2, 3, equal to the number of 
steps of the walker to the left. 

Let us emphasize that our problem can be recast as a purely combinatorial one. 
We assume again for simplicity that v is irrational. Let M n ,k be the number of n-step 
walks starting from the origin, remaining on the left of the wall, and ending at the point 
(n,n — 2k), with the notations Q2.6|) . The numbers M n ,k obey a recursion relation very 
similar to that obeyed by the binomial coefficients (?): 



with the initial condition Mqq = 1, hence in particular Af nn = 1. The recursion ( |3.11| ) 
leads to a double array of integers, shown in Figure 6, that we refer to as a truncated 
Pascal triangle. 

If the slope v is rational, i.e., p c is rational, one has to consider two series of integers, 
M^ki which obey the recursion relations 






(3.11) 




(k > np c ), 
(k < np c ), 
(k > np c ), 
(k < np c ), 



(3.12) 
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again with A/qo = 1- 

As the weight of a finite walk only depends on its endpoint, the integers M n ,k formally 
contain all the relevant information for the computation of the survival probability. More 
powerful techniques allowing to extract this information will be exposed in the following. 



Figure 6: Combinatorial approach: truncated Pascal triangle. 
3.3 Probability flow equation 



The recursion relations (|3.11 ) , ( |3.12| ) for the integers M n k are different from those for 



the binomial coefficients at the crossing edges. This indicates that there is a probability 
flow through the wall. To analyze this probability flow, we define the numbers A k as the 
elements of the truncated Pascal triangle at the beginning of a crossing edge: 



At=K-i*- (3-13) 



Eqs. ( glTD , |U2|) imply that walks are absorbed at the fc-th crossing edge. In 



order to count them with the right probability [see eqs. ( |3.14| ), ( |3.15| ) below], we have to 



put these numbers at the end of the crossing edges, obtaining thus the picture shown in 
Figure 7. In the above example, we get the sequence 1, 2, 5, 19, . . . 

This construction leads to the most important (though elementary) equation of this 



paper. The probability to begin with a left step is p, while A^ walks are lost at the fc-th 



crossing edge, each with a probability p k q n i k , hence the probability flow equation 

Int ± (np c ) 

F ± {n,v)=p- £ A±p k q<-\ (3.14) 

k=l 
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where we have used the equivalence between the inequalities n k < n and k < Int (np c ). 
Taking the n — > oo limit, we obtain the following flow equation 

F ± =p-Y.Atp k q<- k (3.15) 

fe>i 

for the limit survival probabilities F ± , which still depend on p and v. 




Figure 7: The probability flow is characterized by integers A k , living on the crossing 



edges [cf. Figure 5]. In the example, we have A 1 = 1, A\ 



2, At 



5, and so on. 



The flow equations (|3.14 ) and (|3.15 ) will be used several times in the following. 
In particular, eq. (|3.15 ) will be used in sec. 4 in order to derive the second and third 
recursion relations satisfied by the sequences A k . Before we turn to these matters, we 
can obtain some information on the analytic structure of F ± in p. Eq. (|3.15 ) shows that 
the expression 

J2^p k q<- h (3.16) 



k>l 



with q = 1 — p, is a convergent sum of positive numbers for any p in the interval [0, 1]. 
Using the definition ( |3.9|) , we have 



p k q n * - k 



zip) I 9 



k l-FracT(fc/p c ) 



(3.17) 



with 



z(p) = pq 



(3.18) 



As the second factor qi 1-Erac ( fc /p«0 i s bounded, we can concentrate on the first factor 
z(p) k . The maximum of the function z(p) for p in [0, 1] is reached for p = p c . Hence the 
sum (|3.16| ) is a holomorphic function of p in a neighborhood of the origin, and it can be 
expanded term by term as long as \z{p)\ < z(p c ). Suppose now p > p c . Eq. ( |3.17|) implies 
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p k q n k k < p k q™ k for each k, so that 

F ± > p - £ Afpfc^ - * = p - p c . (3.19) 

k>l 

This property gives a proof of the fact, announced before, that F > for p > p c . This 
also shows that the sum ( |3.16|) is not analytic at p — p c , because if it were so the equality 
F ± = 0, which holds for p < p c , could be continued across p = p c . 



4 Arbitrary slopes: recursion relations 

This section is devoted to the investigation of the sequences of integers Af. The slope 
v is an arbitrary number throughout this section, so that we shall omit the superscripts 
±, for brevity. 

4.1 Continuity of the path 

By expressing the continuity of the path of the random walker, we shall obtain the first 
recursion relation defining the integers A^. 

Consider a walk of length that starts to the left and terminates at the end of the 
fc-th crossing edge. There are such walks. Either this walk crosses the wall for the 

first time along the k-th crossing edge (there are A^ such walks), or it has crossed the 
wall for the first time along some earlier crossing edge, say the £-th one. In the second 
case, the full walk is the concatenation of one of the Ai walks that cross for the first time 
along the £-th crossing edge, and of any walk from (rif,£) to (rik,k) (there are 
such walks). We thus obtain our first recursion relation for the sequence of integers A^\ 

Eq. ( |4.1| ) clearly leads to integer values for A^, but not obviously to positive ones. The 
positivity of the A^ can be checked on the first few of them, i.e., 

A 1 = l, A 2 = m-1, A 3 = -(n 1 -l)(2n 2 -n 1 -2), etc. (4.2) 

The recursion relation (|4.1[ ) looks like a convolution. It fails however to be an exact 
convolution, because n^ — ne ^ Uk~£ in general. Nevertheless the difference between these 
integers is small (zero or one in absolute value). 
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4.2 Probability flow 



The second recursion relation for the integers relies on the probability flow equa- 
tion ( 3.15|) , and on the analyticity at small p of the sum (|3.16|) , proved at the end of 



sec. We know that F = for fixed v and small enough p. Hence we can apply 
contour integrals to both sides of the relation 

p= y £A k ^q n "- k . (4.3) 

k>X 

To be more precise, let I > 1 be an integer, and ce(p) be a holomorphic function in 
a neighborhood of the origin, with a Taylor expansion starting as ce(p) = p e+1 + ■ ■ ■ 
Integrating along a small contour around the origin, we have 



T^T? = £ A k I ^j- P k ^~ k = A l + Y,Aj -^ P k ^ k . (4.4) 



Indeed the integral in the two rightmost sides equals for k > £, and 1 for k = £. 

We can take any sequence of functions cg(p), and get a corresponding recursion relation 
which determines the A^. The following three cases will prove useful in the sequel. 

• The most obvious candidate is 

c e (p)=p e+ \ (4.5) 
which leads to our second recursive definition of the A^. 

S k ,i = A k + £ (-l) fc -^"fK ( 4 - 6 ) 



£<k-l, k<n t 



where 5^ is the Kronecker symbol. Again, the integrality of the A^ is obvious, but 
their posit ivity is not. 

Another simple choice reads 

c e (p) =p e+1 q n *- e+1 . (4.7) 

Surprisingly enough, we recover our first recursion formula ( |4.1| ), expressing the 
continuity of the path. This example illustrates the generality of eq. 



Let us come back to eq. ( p.17 ) , which shows that the sum (|3.16|) is almost an entire 



series in z(p), with the definition ( 3-18 ), up to a small positive power of q, with an 



exponent in the range [0, 1]. If we now divide expression ( |3.16| ) by q, we have again 



almost an entire series in z(p), up to a small negative power of q, with an exponent 
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in the range [—1,0]. It is therefore natural to evaluate the contour integrals of 



eq. (|4.4j ) with the weights 



dp dz dp dz 

and — - = —-- (4i 



ci(p) * e+1 c e (p) z*+i(l-p(z)) 

where p(z) is the inverse function of z(p), defined in eq. ( 13.181) . After somewhat 
lengthy computations, which boil down to binomial expansions, we respectively 
obtain 

nk/Pc - 1) _ A ^ T(k/ Pc -n e ) FracU/v))^ 

k\T(k qc / Pc ) ~ Ak t 1 ^-i)\nk qc /p c + l-n l + ^ 1 ftaC ™ ) ^ 

mPC+1) - * + .. Frac^A, (4.9) 



k\ T(kq c /p c + 2) fri (k-£)\T(kq c /p c + 2-n i + 



The recursion relations (|4.9| ) seem rather complicated, so that even the integrality of 
the Ak is not obvious. Eqs. ( f4.9| ) are nevertheless of interest. Indeed all the coefficients 
in the sums have the same sign, since the are positive. We thus obtain the following 
bounds 

r(k/ Pc - 1) r(fc/ Pc + i) 

k\T{kq c /p c ) ~ k ~ k\T(kq c /p c + 2) 1 ' } 

The case where p c = 1/N, with N an integer, is of special interest. In this situation, 
ijpc = £N is an integer for any £. Re-introducing the ± superscripts for a while, we 
thus have Frac + (£/p c ) = 1 if we are interested in F~, and Frac~(£/p c ) = if we are 
interested in F + . As a consequence, the A^ saturate the lower bound of eq. ( |4.10|) , while 
the A£ saturate the upper bound. We shall recover these properties in sec. [0| with more 
powerful techniques. 

In the general case, we set 

A * = B * fcl T(kq!/ P pc + 1) ' ^ 



Eq. ( [4.1 0| ) shows that the prefactors B k are bounded. In other words, using again 



Stirling's formula, we obtain the estimates 



A±k * B * (2^)1/2 fa (* » 1)- (4-12) 

The prefactors obey the simple bounds 

q c < Bt < (4.13) 
Qc 

We shall come back in sec. ^| to these prefactors, which are non-trivial in general. 
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4.3 Duality 



The recursion relations derived in the previous section are based on the fact that F = for 
small enough p. We now aim at using the vanishing of F on the full interval < p < p c . 
We call the following approach duality, as in the beginning of sec. 0, since it relates 
the survival probability in presence of the obstacle with slope v to that with slope —v. 
In sec. [5.3| we shall give a simple combinatorial argument explaining why the slopes v 
and —v, i.e., p c and q c , are related, at least for rational v. 

We start again from the identity (|4.3|), 



p-$> fc pV*- fc = (0<p< Pc ), (4.14) 
fc>i 

and the similar equation for the slope —v, with integers Ak, 

p-J2 A k p k q^- k = (0 < p < q c ). (4.15) 

fc>i 

By substituting q = 1 — p for p in this last equation, we get 

q - £ A k q k p^- k = (p c <p<l). (4.16) 

k>l 

Multiplying eqs. (|4.14|) , ( |4.16|) , we obtain 



p-J2 Mp k q nk ~ k )[q ~ E A k q k p nk ~ k = (0 < p < 1). (4.17) 

fc>l / \ k>l J 

For the rest of this section, we assume that v is irrational, while the case where v is 



rational is treated in sec. |5.3j . For v irrational and k > 1, the integer n k corresponding to 
the fc-th crossing with the wall at slope v is defined by the inequalities < p c n k — k<p c , 
while the integer n k corresponding to the k-th crossing with the wall at slope — v is 
defined by the inequalities — q c < p c n k — (n k — k) < 0. Hence the fractional parts of p c n k 
and p c n k live in the disjoint intervals [0,p c ] and [— q c , 0]. Consequently, the set of the n k 
and the set of the n k have no integer in common. Note that n k and n k are both larger 
than 1. 

Conversely, for n > 1, there is a unique integer k n such that — q c < np c — k n < p c . We 
set 9 n = np c — k n . Equivalently, with the definitions (|2~8|), we have 



k n = l + Int((n-l> c ), 9 n = -g c + Frac((n - l)p c ). (4.18) 

Let J + be the set of integers n > 1 such that 6 n > 0. For such an n, there is an 
integer k > 1 such that n = n k , namely k = k n , and we have 

nel + : Q n = p c (l — Frac (~J J ■ ( 4 - 19 ) 
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Similarly, let I be the set of integers n > 1 such that 9 n < 0. For such an n, there is an 
integer k > 1 such that n — rik, namely k = n — k n , and we have 

?7i kri \ \ 



n G I : = — 9c 1 — Frac 



(4.20) 



The above two sets are a partition of the integers n > 1. We can therefore define a 
sequence of integers A n for n > 1 by 

f A fcn if n e /+, 
[A n _ fcn if nel~. 



(4.21) 



A last geometrical observation is in order. It follows from the definition of the sequence 
k n that k n < k n+ i < k n + 1. Hence the sequence 



(4.22) 



defines a canonical walk, shown in Figure 8. 

This walk is closest to the wall in a well-defined sense: the point (n, n — 2k n ) is the 
endpoint of an edge crossing the wall, either from left to right or from right to left. 
Equivalently, the next step of the canonical walk is to the left if it starts to the right of 
the wall, and vice versa. We identify the symbol X n with the polynomial 

X n =p k "q n - k ", (4.23) 

whose total degree is n, and whose difference of partial degrees in p and q gives the 
position X n , as in eq. (|4.22|) . 




Figure 8: The canonical walk. 

With the notations ( [OTP , QE2g) , eq. (gT7|) can be recast as 

l-q- J2 A n XA lq-J2 A n X n ] =0 (0<p< 1) 
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(4.24) 



Let us expand the product. It is easy to check that if n and m have different /-signs, 
then X n X m = X n+m . Moreover, if n £ I + , pX n = X n+ i, and if n £ I", qX n = X n+ \. 
Finally, pq = X 2 . We set A\ — 1, and Iq = I + U {1} and J^" = I U {1}. We obtain after 
some algebra 

£3 n X n = ~A m A n X n+m (0<p<l). (4.25) 



n>l 



A heuristic term-by-term identification leads to a third recursive definition of the integers 
A n '. 

A e = J2 (£>2). (4.26) 



m+n=l 



An algebraic proof of this relation will be given in sec. 5.3 



5 Rational slopes: algebraic approach 

In this section we investigate in further detail the outcomes of the techniques exposed 
in sec. f|, in the case where v is rational. Algebraic functions play a central role in this 
analysis. 

We write a rational slope in terms of two relatively prime integers M, M > 1 as 

~ M-M MM , . 

N = M + M, v= - , Pc = j^, ^=jv- ( 51 ) 



The integers n^, defined in eq. ( |3.10|) , are skew periodic, in the sense that 

nt+M = nt + N, (5.2) 

and we have 



n. 



n+ (k = l,...,M-l), n M = N } n + M = N+l. (5.3) 



We recall that, for each rational slope v, we have to consider two different sequences of 
integers, Af, associated with the survival probabilities F ± . 

5.1 Continuity of the path 



As a consequence of the property (5.2) of the integers n k , the recursion relation ( 4.1 



can be split into M blocks, one for each value of k modulo M. 
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We introduce the following formal power series in t: 

Ft(t)=Y. A± KM +k t K (* = 1,...,M), (5.4) 



K>0 



<%W=X( K fjijV (* = 1,...,M), (5.5) 



,. „ v KM + k-l 



K>0 
KM+k-£>0 



Eq. ( |4.1| ) becomes a set of M coupled linear equations, of the form 

M 

E GkfiWV) = Gf (*) (A: = 1, • • • , M). (5.7) 

The solution of this linear system yields the functions F fc (i), which are generating func- 
tions for the integers . 

So far t is a formal expansion variable. There is however a natural choice for it, 
namely 

t = p M q M = p M (l - p) N ~ M = z{p) M , (5.8) 

so that the survival probabilities F ± can be expressed in terms of the functions F^{t). 
Indeed, with the definitions (|5.8j ) of the variable t and (|5.4j) of the functions Fj^(t), the 
probability flow equation ( |3.15| ) can be recast as 

M 

F ± =p _Y / p k q n *- k F±(t). (5.9) 

k=l 

We now show that the functions Gf(t) and G^ e (t) are algebraic functions of t. These 
functions are special cases of the functions 

««w= («») 

for i,j non-negative integers. Our starting point is the following contour integral repre- 
sentation of the binomial coefficient: 

fl\ f dz (1 + z) 1 r du 1 



J I I 2vri z J+1 



/ 2vri u J +m-uY- J + v (5 ' U ' 



where the integral is along a small contour around the origin. The second expression 
is obtained by means of the change of variable z = u/(l — u). Applying eq. ( |5.11| ) to 
(/, J) = (KN + j, KM + i), and summing over K, we get 

2?ri u M(l - U )N-M _ t 



H dt) = f ., M( \_Jn- M 7 ■ ( 5 - 12 ) 
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Now, we need to know which roots of the denominator are inside the integration 



contour. The geometric sum over K is convergent for \t\ < 
other hand, the polynomial equation 



u 



Mi 



u 



\N—M 



U M (1 _ u) N-M = f 



has M 'small' roots u n , with a — 1, ... , M, such that 



t l/M u a-l 



(t -> 0), 



with 



The small roots can be continued in t. For t real in the range < t < t c , with 

_ M M M {N-M) N ~ M 



N N 



On the 



(5.13) 



(5.14) 



(5.15) 



(5.16) 



the M small roots remain inside the circle \u\ = p c , while the other N — M roots of 
eq. (|5.13 ), the 'large' ones, lie inside the circle |1 — u\ = q c . The two circles in the it- 
plane are tangent to each other at u — p c . For p < p c , one of the small roots is p itself, 
while for p > p c , one of the large roots is p itself. Figure 9 shows a plot of the roots 
of eq. (|5TT3| ), for M = 4, iV = 11, and p = 1/2 > p c = 4/11, i.e., t = 1/2048, while 
t c = 108/3125. The root p = 1/2 is shown as a large empty symbol. 



0.E 



0.4 



-0.4 




-0.E 



-0.6 



0.6 

Re u 



1.2 



Figure 9: Roots of the polynomial equation (|5.13|) for M — 4, N — 11, and p = 1/2. 
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For \t\ small enough, the contour of eq. (|5.12 ) contains exactly the M small roots u a , 
so that we have 

M i 

Each of the is an algebraic function of degree iV, since it obeys eq. ( |5.13|) . As a 

consequence, Hij(t) is an algebraic function of degree (J^), for any non-negative integers 
i, j. Indeed, choosing one branch of the function Hij(t) amounts to choosing M branches 
of u a (t) among N. 

The algebraic function Hij(t) is singular when two roots of the polynomial equa- 
tion ( f5.13| ) coincide, in such a way that the contour of eq. (|5.12| ) gets pinched. The only 
non-trivial singularity corresponds to t = t c , where the real positive small root u = u\ 
merges with the root u = p into a double root at u = p c . In the vicinity of this point, we 
introduce the notation 

p = Pc + 5p, t^t c (l--^(Sp) 2 ) . (5.18) 

For 5p > 0, the first root Ui has a correction in Sp, while the other M — 1 roots are 
regular in t, so that their corrections are of order (Sp) 2 : 

Ut^pc-Sp, u a = (u a ) c + 0((Sp) 2 ) (a = 2,...,M). (5.19) 

We thus have M — Nui ~ iV Sp, hence the estimate 

1 1 / t c \ 1/2 

Hii ® ~ plqi' l NSp ~ p^ 1 \2Np c q c (t c -t)) ' (5 ' 20) 

This inverse-square- root singularity can alternatively be obtained by means of eq. ( 2.20j) , 
from the large- K behavior of the coefficients of the series (|5.10|) . 

We have thus shown that the functions G^(t) and G^ e (t) which enter the linear 
equations ( pT7|) are algebraic in t, with degree (J^j- They possess branch points of the 
form (|5.20j) at the critical point (t = t c ) . Unfortunately, these properties are not sufficient 
to evaluate the critical behavior of the functions Fj^(t) near t c . Indeed, it can be checked 
that the linear system (|5.7| ) becomes very singular at t — t c : the most singular part of 
the matrix G^ £ (t) as t — > t c is a matrix of rank one. 

It is worth noting that the small roots of the polynomial equation ( 5.13Q also play an 
important role in recent works [f§ devoted to large deviations in exclusion models. 



5.2 Probability flow and algebraic trick 

An algebraic treatment can also be applied to eq. ( |5.9| ). Consider a fixed p > p c , so 
that t < t c . If in eq. (|5.9|) we replace p by any of the small roots u a of eq. (|5.13| ), the 
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r.h.s. vanishes. Indeed this expression has been shown below eq. ([3.18|) to vanish for 
\z(p)\ < z(p c ), and we have indeed |^(w a )| = \ z (p)\ < z ipc)- Hence 



M 

+ 



u a = J2 " ^T k (a = 1, . . . , M). (5.21) 

k=l 

This linear system of M equations determines the Fj^(t). Eq. ( |5.9|) then leads to the 
survival probabilities F ± . 

This algebraic trick is quite powerful for numerical purposes. Indeed, for given p, 
it only involves solving the polynomial equation (|5.13|) and the linear system (|5.21| ). 
Figure 10 shows a plot of the survival probabilities F against the slope v, for different 
values of p, indicated on the curves. For given p, the algebraic trick has been applied 
numerically for all rational slopes having a denominator iV < 40. The numerical values 
of F^ 1 are thus obtained. For a rational slope v , F + (respectively, F~) can be read as 
the ordinate just after (respectively, just before) the discontinuity of the curve F(v). 
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Figure 10: Plot of the survival probabilities F^ 1 against v, for several values of p, 
indicated on the curves. The apparent v — > 1 limit of F(v) is F~(l) = p [see eq. (|7.5|) 1, 
whereas the jump at v — 1 to F + (l) = 1 is not visible. 

The algebraic trick also demonstrates that the Fj^ are generically algebraic func- 
tions of degree deg t (M, N) in t, while F ± are generically algebraic functions of degree 
deg p (M, N) in p, with 



deg t (M,N) 



(5.22) 
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Elimination theory can be used to write down explicitly the algebraic relations between 
p or t and the F k or F ± . Various examples will be given in sec. 0. 

If we are interested only in F ± , there is an alternative routed]. Instead of computing 
the F k , we can solve a transposed linear system, whose meaning is the following. Suppose 
that p is replaced by one of the small roots u a to count the weight of the walks. Then 
the crossing edges that end on the wall get the correct weights, while the other ones are 
slightly wrong, because nf/k is only approximately equal to N/M. By taking a linear 
combination of the weights corresponding to the u a , we can give the right weight to any 
crossing edge. Indeed, given the u a , we can determine numbers such that we have, 
for any k > 1 

M ± ± 

Y, Dtu k a (l - u a ) n k ~ k = p k q n k - k . (5.23) 

This is possible because eq. ( |5.23| ) has to be imposed only for k — 1, ... , M, as a conse- 
quence of eq. (|5.2|) . We can substitute eq. (|5.23| ) into eq. (|3.15| ), and use again the fact 
that the r.h.s. of eq. ( |3.15| ) vanishes when we substitute u a for p, obtaining thus 

M 

F ± =p-Y D±u a . (5.24) 

a=l 

We close up this section by investigating the relationship between F + and F~ . Using 
eq. flOl) , we can recast eq. (|5.9| ) as 

M-l 
k=l 

M-l 

F- = p -Yp k q n *- k F k -tF M . (5.25) 
fc=i 

Since the above equations determine the F k , a comparison between both equations yields 
or equivalently 

F + = T^Tj^> ^ + = T J T^' (k = l,...,M). (5.27) 
Finally, eqs. ( |3.8| ) and ( |3.7[ ) respectively become 

Q(v) = tF M (5.28) 



and 



= = = ^mF + = tF+F~. (5.29) 



1 This approach was suggested to us by V. Lafforgue. 
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5.3 Duality 



We now turn to an investigation of the relationship between the survival probability 
associated with the rational slopes v and —v, by exploring the duality approach of sec. [O] 
in the rational case. We notice that substituting v for — v amounts to exchanging the 
roles of the integers M and M in the definitions (|5.1|) . 

We start with the following combinatorial observation. Let (t = n, x = n — 2k) be a 
point with integer coordinates on the wall. The number of finite walks starting at the 
origin, remaining on the same side of the wall and ending at (t, x) is the same for left 
walks and for right walks. There is indeed an obvious one to one correspondence between 
both sets of walks, depicted in Figure 11. The second walk (contributing to right walks) 
is obtained from the first one (contributing to left walks) by a rotation of angle 7r around 
the mid-point (t/2,x/2). 



Figure 11: Combinatorial proof of the identity ( |5.30|) . 

This observation implies the following identity between generating functions 



(5.30) 



In particular, the probability Q that a walk crosses the wall and does so for the first 
time at a point with integer coordinates, as given by eq. ( |5.28| ), is the same for right and 
left walks. We shall not use this information for a while, but rather derive many other 
identities, including eq. ( |5.30| ), by employing only algebraic means. 

Following the same strategy as in the irrational case, we start from eqs. fl4.14| ), ( |4.15| ), 
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which now read 

M 



p-EpV r ^W= (0<p< Pc ), (5.31) 
fc=l 

M 

p-E P V^(i)=0 (0<p<g c ), (5.32) 



fc=i 



with t = p M q M . If we now substitute q = 1 — p for p in eq. ( |5.32| ) (this turns t into t), 
and then multiply it by eq. ( |5.31| ), we get 



/ M \ M ~ ~ \ 

(l-q-J2p k q ni - k F±j U-E^-^j =0 (0<p<l). (5.33) 

We again use the canonical walk, but now only a finite piece of it. For definiteness, 
we concentrate our attention on the FjT, for both slopes v and —v. We first observe that 
n M = ^ = iV + 1, so that, if we define 

F + = 1 + tF+, F + = l + tF±, (5.34) 

eq. ( |5.33| ) becomes 

/ M -! \ / A^- 1 _ \ 

U-qF+- Y.P k( f +k ~ kF k) l-^o + - =0 (0<p<l). 

(5.35) 

We now change notation for the unknowns inside the sums, just as in the irrational 
case, labeling the F + functions by the value of or n^, rather than with k. We employ 
the notations ( |4.18| )-( |4\23| ), albeit with the sets J 1 * 1 now being a partition of the integers 
£ = 2, . . . , N — 1. Eq. Q05D reads 

1 " <lK - E ) \ ±-pK- E ] = (0 < V < I)- (5.36) 

The strategy for expanding this product is the following. We substitute t for p M {l — 
p) as often as we can. For fixed t, we thus end up with a polynomial in p of degree 
less than N. Such a polynomial has to vanish identically, because it has iV roots. Indeed, 
the first factor in eq. ( |5.35[ ) vanishes when p is equal to any of the M small roots u a of 
the polynomial equation (|5.13 ), while the second factor vanishes on the other M roots. 



Now consider I e J + and f G I~ . If £ + £' < N, then X e X e = X e+t . If £ + £' = N, 
then X e X e = t. If £ + £' > N + 1, then X t X e = tX w ,_ N . The equality £ + £' = N + 1 
never occurs. We also have pXj> = X^ +1 if £ + 1 < N and pX N -\ — t if iV — 1 G J + , 
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and symmetrically qXpj = Xgi + \ ii £' + 1 < N and qN-i = t if N — 1 G I~ . Finally, 
pq = p(l — p) = X%- All these properties can be checked by combining the inequalities 
defining kg and kg> [see eq. ( [4.18| )1. All the Xgs have different degrees ranging from 2 



to N — 1, so that, together with any two of the three polynomials 1, p and 1 — p, they 
form a basis of the polynomials of degree less than N. We concentrate for a while on 
terms of degree less than two in p. We are left with the sum of a function of t and of 
— (1 —p)Fq — pFq , which has to vanish identically This implies in particular Fq~ = F + , 
which is just the identity ( |5.30|) , that we have proved in the beginning of this section by 



combinatorial means. We set 

Ft = F+ = F+, (5.37) 

and Iq = I + U {1} and Iq = I~ U {1}. Just as in the irrational case, we meet many 
simplifications, and we are left with 

AT-l 

1 - F\ ~ E Ftx e + Yl FtFtt p ^X w ,_ Np{e/l) = 0, (5.38) 

with the convention X = 1, and with 

,„ „, s ( if i + £' < N, . . 

p(£,t) = { , 5.39 

^ V ' ; \ 1 if £ + £' > N. K ' 

A term-by-term identification then leads to 

Ft = l + t ]T Ft Ft (5.40) 

e+e'=N 



and 

Ft= E FtFU^ (£ — 2, . . . , N — 1). (5.41) 

£'+£"=£ mod N 

Eqs. (|5.40|) and (|5.41| ) provide N — 1 equations for the N — 1 unknowns F~t ■ When 
t goes to 0, eq. ( |5.40| ) becomes trivial, while the other ones become identical to the 
recursion relation ( |4. 26|) . As any irrational slope v can be approximated by a rational 
one that leads to the same coefficients Ai up to an arbitrary given £, we thus obtain a 
proof of the third recursion relation ( |4.2(j| ) by algebraic means, without any recourse to 
analysis. 

6 Critical behavior 

This section is devoted to the analysis of the 'critical behavior' as p — > p c of the survival 
probabilities F ± and of related quantities in the convergent regime (p > p c ), for a fixed 
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slope v, either rational or not. This investigation also yields quantitative predictions 
concerning the amplitudes C ± (v) in the marginal regime (p = p c ) and the prefactors 
b ± (np c ) in the large-deviation regime (p < p c ). 



6.1 Convergent regime 

We consider first the convergent regime, for a rational slope. The behavior ( 4.12 ) of 
the integers implies that the generating series Fj^(t) have square- root singularities at 
t = t c , as shown by eq. (|2.20|) . We again employ the notations (|5.18|) , and we define the 
positive amplitudes $^ such that 

Ftit) ~ (if )c - H \Sp\ (k — 1, . . . , M, 5p —> 0), (6.1) 

with (Fj^) c = Fj^(t c ). Similarly, the survival probabilities F ± are expected to vanish 
linearly as p — > p c + 0, with amplitudes $ ± : 

F ± &$ ± 5p (5p^+0). (6.2) 

The critical amplitudes $ ± and $^ can be determined by means of the algebraic trick 
as follows. By differentiating eq. (|5.9|) , we obtain 



^ =1 -£^ {{p-^r) Ft+ ^\- (6 - 3) 



Now, by inserting the limits 



dF± *± dF k 

0, ^ -<IV (p -/v o). 



(6.4) 



dp dp 

into eq. (6.3 ), and using eq. (|3.10 ), we obtain two different expressions for the amplitudes 



M + \ M 

n,-k-l 



$± = 21 + EPc9c fe (l - Frac^/j9 c ) (F ± ) c = 2$> C V (6.5) 



k=l / k=l 



The middle side of eq. ( |6.5| ) gives $ in terms of the (Fjr) c , which can be obtained 
by the algebraic trick at p = p c . The rightmost side of eq. ( |6.5| ) then allows to determine 
the $± by means of a second use of the algebraic trick. Indeed the behavior ( 5.19 ) of the 
roots u a implies du a /dp — > —S a ,i in the p — > p c + limit. As a consequence, the ^ can 
be determined from the following linear system of M equations 

s: m 

-f^ = EWc(l - Wc)"' 4 ^ (« = 1, . . . , M). (6.6) 
z fe=i 
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It turns out that the amplitudes $ ± and govern the behavior of various quantities 
for a rational slope. First, eq. ( [2.20D yields the asymptotic behavior of the integers A^ 
for large k. We thus recover eq. ( f4.12| ), and we obtain expressions for the amplitudes B^, 
which asymptotically only depend on k modulo M: 

B^M+k « Np c q c (t c ) h / M $± (A; = 1, . . . , M, K > 1). (6.7) 

This result can be recast as follows. Since p c = M/N is rational, the fractional part 
Frac ± (fc/p c ) = Frac ± (/cA^/M) takes M rational values with denominator M, and only 
depends on k modulo M. Conversely, any function of k modulo M can be considered as 
a function of Frac ± (/c/p c ). Hence eq. (16.71) is equivalent to 

Bt « B ± (k/p c ), (6.8) 



i.e., 



where B ± (x) are periodic functions, with unit period. Since the functional form of 
eq. (|6.9|) holds independently of the slope v, provided v is rational, we make the reasonable 
hypothesis that the law (|6.9|) also holds when v is irrational, albeit with one single periodic 
function B(x). 

The bilateral sequence of integers A^, introduced in eq. (|4.21|) , can be argued to 
exhibit an asymptotic behavior analogous to eq. (|6.9|) : 

^^ w i^F ^ r (n>>1) ' (6 ' 10) 

where the arguments 8 n of the periodic functions /3 ± (^) have been defined in eq. (|4.18| ). 
The estimate (|6.10 ), with its periodic functions (3 ± {9) 1 involves p c and q c in a symmetric 
fashion, contrary to eq. (|6.9|) . Furthermore it encodes simultaneously the asymptotic 
behavior of the sequences associated with both slopes ±v. These are advantages of 
the duality approach. Eq. fl4.19|) also implies the following relationship 

B ± (x) = (p c qf l^)~ e (3 ± {6), 9 = p c (l -x) (0 < 9 < p c ), (6.11) 

between the periodic functions of the estimates ( |6.9|) and ( |6.10|) , with corresponding 
arguments. Figures 12 and 13 respectively show plots of the periodic functions B(x) and 
(3(9), again for the golden slope given in eq. ( |2.23| ). 



6.2 Marginal and large-deviation regimes 

Quantities pertaining to the marginal regime (p = p c ) and to the large-deviation regime 
(p < p c ) can also be characterized in terms of the critical amplitudes $ ± and for a 
rational slope. 
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In these two regimes, the limit survival probabilities F vanish, so that the probability 



flow equation (|3.14 ) becomes 

F ± {n,v) 



fc>Int ± (np c )+l 



(6.12) 



1.2 



1.1 



X 

PQ 



1.0 



0.9 



0.8 




0.2 0.4 0.6 

X 



0.8 



1.0 



Figure 12: Plot of the periodic amplitude B(x), for the golden slope. 
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Figure 13: Plot of the periodic amplitude (3(9), for the golden slope. 
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In the marginal regime (p = p c ), the behavior of F ± (n,v) for large n can be derived 
by inserting the estimates (|4.12| ) and (|6.7|) into eq. (|6.12|) . The terms in the r.h.s. of 
that equation decay as the power law k~ 3 / 2 , so that the result is rather insensitive to 
the details of the sequences A k . Indeed the amplitudes $^ only appear through the 
combinations which enter the rightmost side of eq. ( |6.5|) . We thus recover the general 
result ( |2.29| ), with an amplitude 

/1 „,2\l/2 



C {v) 



/'<■</■ V /2 $± 



(6.13) 



in agreement with eq. ( |2.35| ) 




V 



Figure 14: Logarithmic plot of the critical amplitude C(v) against v. 

The algebraic trick at the critical point thus allows a numerical determination of the 
amplitude C{v) of the law (|2.29|) . Figure 14 shows a logarithmic plot of this amplitude, 
obtained from data at rational slopes, along the lines of Figure 10. The curve is centrally 
symmetric, as a consequence of the duality symmetry ( |2.30| ). Figure 15 shows a similar, 
albeit more appealing, plot of the combination 



y ± {v) = lnC^) + ^ In — = lnC^u 



1, l+v 

- in . 

2 1 - v 



2 p c 

This quantity obeys the following bounds, shown as dashed lines in Figure 15: 



(6.14) 



In 2 , , , In 2 , . 

< y (v) <— In(l-u) 



2 - a v ' ~ 2 

x In 2 , , N In 2 
ln(l + u) < y ± (v) < 



(-1 < v < 0), 
(0 < v < 1). 



(6.15) 
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The results ( [7. 12| ) and ( 7.27 ), to be derived later, show that the bounds ( p,15| ) are sat- 
urated by the slopes with either M = 1 or M — 1. These bounds lead to the limit 
values 

/ \ In 2 
y(±l) = ± — . (6.16) 
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Figure 15: Plot of y(v), defined in eq. (|6.14j) , against v. 



The investigation of the large-deviation regime (p < p c ) is more involved. The sum 
in the r.h.s. of eq. ( |6.12 ) is now exponentially convergent, so that it is dominated by the 
fine structure of the near the lower bound, i.e., for k rs np c . We recover after some 
manipulations the general result Q2.21|) , with amplitudes which only depend on n modulo 
N: 



u KN+n ~ 



Nq c 



M 



^ g l-FracT(fc/ Pc ) e [(n-k/ Pc )-NB(n-k/p c )]S(v) ^Qk/M^±^ (6.17) 



for n = 1,...,N and K ^> 1, where is the Heaviside step function, and S(v) the 
entropy function of eq. ( |2.13| ). This result agrees with the existence of a periodic prefactor 
to the law ( |2.21| ), since any function of n modulo can be viewed as a function of 

Frac ± (np c ) = Frac ± (nM/A^). 



7 Further results in specific cases 
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7.1 The limit case v = 1, i.e., p c = 



In this limit situation, quantities with the superscript + are equal to the same quantities 
for the whole unconstrained set of random walks. Quantities with the superscript — are 
only constrained by the condition that the first step is to the left (e\ = —1). 

The general formalism of this paper agrees with these considerations. Indeed, any 
value of p > corresponds to the convergent regime. Eq. (|2.7|) implies that the one-time 
distribution functions read 

P + {n,v) = l, p-(n,v) = l-q n (n>l), (7.1) 



in agreement with the value S = — Inq of the entropy function [see eq. ( p.l3|) ]. We have 
therefore 

f + (z,v) = -^—, = Vz^' ( 7 - 2 ) 

l Z X z 

so that 

F + (n,v) = l, F~{n,v)=p (n > 1). (7.3) 
The integers M^k are a ^ so remarkably simple: 

^ = (l-i) ( n > i ^ k > 1 )- ( 7 - 4 ) 

The limit survival probabilities read therefore 

F+ = l, F~=p, n = g, (7.5) 

which can be recovered by elementary considerations: we have II = Prob{F = 1} = 
Probjei = 1} = q and F + = 1, hence F~ = F + — 11 = 1 — q = p. 

7.2 The limit case v = — 1, i.e., p c = l 

This situation is the dual of the previous one. Quantities with the superscript — vanish 
identically, while quantities with the superscript + have the single contribution of the 
walk with all steps to the left (e n = — 1 for all n). 

Any value of p < 1 now corresponds to the large-deviation regime. Eq. (|2.7|) leads to 



P+(n,v) =p n , P~(n,v) = (n>l), (7.6) 
in agreement with the value S — — hip of the entropy function. We have therefore 

f + (z,v) = -^—, /-(*,*) = 0, (7.7) 
1 — pz 
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so that the survival probabilities read 

F + (n,v) =p n , F-(n,v) = (n > 0). (7.8) 
The only non-zero integers M^k are 

Af+ n = 1 (n > 0). (7.9) 

7.3 Rational slopes with M = 1, i.e., v = 1 - 2/N or p c = 1/iV 

The class of rational slopes with M — 1 owes its simplicity to the fact that only one 
unknown function, F^(t), is involved in the algebraic trick approach. 

For M = 1 and iV > 2, we have v = 1 - 2/N, p c = 1/N, q c = (N - 1)/N, t = 
pq N_1 , and t c = (N — l)^ -1 /N N . In the convergent regime (p > p c ), the polynomial 
equation (|5.13|) has only one small root, u\, which is in the interval < U\ < p c . Eq. ([5.21D 
has the explicit solutions 

*T = ^, Ft = -* . (7.10) 

t t(l — U\) 

We thus obtain 

F =p-ui, F + = - , Q = ui, n = — ■ . (7.11) 

1 — Ml 1 — Mi 

We notice that ui(t) and F^(t) are algebraic functions of degree iV in t, while Mi(p) and 
F ± {p) are algebraic functions of degree N — 1 in p, in agreement with eq. ( |5.22|) . 

The critical behavior of the quantities given in eqs. ( 7.10|) , ( 7.11 ) is characterized by 
the amplitudes 

/ N \ N ~ 1 ( N xiV 



^1 / Ar 1\)V-1 ' ^1 



(JV-l)JV-i' a (jv_i)jv+i' 
The integers can be obtained in several ways. The easiest route consists in eval- 



uating them as contour integrals of the generating series F{ (t). The expressions ( [7.101 ) 
yield 

/ dt [dm 1 - Nm 

2%iu k (1 — Mi 



Ak J 2mt k+1 Ul t 2mu k (1 - u,\ Nk ~ k ^ 



+ • dt u x r dug 1 - Nui 

k ' f 2nit k + 1 l-u 1 J 2mu k {1 - Ul ) Nk ~ k + 2 ' 1 J 
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The integrals over t are transformed to integrals over U\ by means of the identity t = 
ui(l — ui)^ -1 , whereas the latter are linear combinations of the integrals over u which 
enter eq. flS.lip . We thus obtain the integers as ratios of factorials: 

A - ( ivfc - 2 ) ! A+ ( Nk y- (7U ) 

k k\(Nk-k-l)V k k\(Nk-k + l)V y ' 

These formulas can also be obtained combinatorially, by means of the Raney principle 0. 
They agree with the observation, made below eq. (|4.10|) , that the saturate the lower 
bound, while the saturate the upper bound. The corresponding amplitudes Bf read 
asymptotically 

N — 1 1 N 

Bt«- =57rT . (7.15) 

The situation of an immobile obstacle, i.e., of a vertical wall, with slope v = 0, 
corresponds to M — M — 1 and N = 2. In this case, we have p c = 1/2 and U\ — 1 — p 
for p > 1/2, so that the results (f7T0|) -(pID can be further simplified to 

Fr = F? = ± F~ = 2p-l, F + = ^, 

p p z p 

Q = 1-p, n-tfr- 1 " 1 -"). (7.16) 

p 

The discontinuity II is maximal for p = l/y/2, where it equals 

n max = (V2- l) 2 = 0.171573. (7.17) 

This value is the largest possible discontinuity met in this problem, except the trivial one 
at v = 1, namely II = q [see eq. (|7.5|) 1- 

This example can be alternatively investigated ab initio. Indeed, the integers J\f n k 
which obey eq. ( |3.12j ) can be written as differences of binomial coefficients: 

Eq. ( |3.13|) then yields, in agreement with eq. ( [7.14] ) for iV = 2, 

(2k-2\ (2k-2\_ {2k-2)\ 
Ak ~ ^^-[k-l)~{ k ) ~ k\ (jfc-l)! -Cfc - 1 ' 

where the are the Catalan numbers fl. 
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7.4 Rational slopes with M = 1, i.e., v = -1 + 2/N or p c = 1 - l/N 

The class of walls with M — 1, i.e., a slope v = 2/N — 1, gives a nice illustration of the 
duality approach, exposed in sec. |5\3| . 

In this situation, we have M = N — 1, p c = 1 — l/N, q c = l/N, t = p N ~ 1 q, and 
t c = (N — l)-^ -1 /N N . According to the general rules of sec. |5l| we introduce N — 1 
functions F±, . . . , F N _ V The function F~^ pertains to both v and —v, while all the other 
ones only pertain to v. We have Iq = {1} and 1$ = {1, . . . , N — 1}, so that the duality 
equations (|5.40|) , (|5.41|) respectively read 

Ft = l + tFp%-i, 

Ft = F^fU (k = 2,...,N-l). (7.20) 

1_ / l\ k 

So, we have F k = IF 1 J for k — 1, . . . , N — 1, and 

Ft = l+t(Tt) N . (7.21) 

The solution Ft(t) is the branch of eq. ( |7.21 ) which is regular at t = 0. A comparison 
of eq. ( |7.21| ) with the identity t = p N ~ l q shows that Ft = 1/p. We thus have 

F + k =- k (k=l,...,N-l). (7.22) 

In order to obtain a complete combinatorial description of the problem, we must recall 
the relation between the functions Ft and the functions F k . We have 

n+ = k + l (k = l,...,N -2), n+_ 1 =N + l, (7.23) 

so that F~t = 1 + tF+ = 1 + tF? and F~t = F fc + _ 1 for k = 2, . . . , N - I. We thus obtain 
the remarkably simple result 

= (k = l,...,N-l). (7.24) 
Eqs. Q and (|06| ) - (|09|) then yield 



F k =- (k = l,...,N-l) (7.25) 

p K 



and 



F- = N P - {N -l), ^ = ^-(^-D , Q = ,, n= ^p-(iV-l)) (726) 

Eqs. (|7^)-(|7^B|) show that the F±(t) are algebraic functions of degree iV in t, just 
as p(t), while F ± (p) are rational functions, i.e., algebraic functions of degree 1 in p, in 
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agreement with eq. (|5.22|) . Finally, the results ( 7.16| ) can be recovered by setting N = 2 
in eqs. 

The simplicity of the results ( |7.24| )-( ]7.26| ), especially for F~, suggests that they can 
be obtained by elementary means. This is indeed the case Q Fix p > p c = M/N. The 
probability to start to the right and to touch the wall is q = 1 — p, because if the first 
step is to the right (probability q), then the walk crosses the wall with probability one. 
Let R be the probability to start to the left and touch the wall. Then, by definition, 
F~ = p — R. The main point is as follows: a walk touching the wall touches it with 
probability one at least once at a point with integer coordinates. If the walk starts to 
the right this is clear, because M — 1, so passing through an integer point is the only 
way to touch the wall coming from the right. A walk starting to the left and touching 
the wall either does it for the first time at an integer point (and we are done), or it 
crosses it at some point from left to right, and then with probability one touches the 
wall again later coming from the right, this time at an integer point, because this is the 
only way to touch from the right. In particular, R is also the probability to start to the 
left and touch the wall at an integer point. Now, if we pick an integer point on the wall 
and compute the probability that a walk starts to the left and passes though this point 
versus the probability that a walk starts on the right and passes though this point, the 
p-dependent weight is the same for both, and what remains in the ratio is the quotient 
of two binomial coefficients, giving M. Although a walk can pass through several integer 
points, as this ratio is the same for any integer point, this is enough to show that R is 
exactly M times the probability to start to the right and cross the wall. The latter is q, 
so R = Mq and F~ = p — R = p — Mq = Np — M, in agreement with eq. (|7.26| ) . 

To finish the combinatorial analysis, we show that Q = q, again in agreement with 
eq. ( |7.26|) . This is a consequence of an observation made at the beginning of sec. [5]3 



and summarized in eq. ( [5.30| ). The probability Q that a walk crosses the wall and does it 
for the first time at a point with integer coordinates is the same for right and left walks. 
But for right walks Q = q, because any walk starting to the right has to cross the wall, 
and can only do so at an integer point. 



The critical behavior of the quantities given in eqs. (|7.24| )-( |7.26| ) is characterized by 
the amplitudes 

N k+i / j\r s k+2 



N 2 /N-l\ 1/2 N 



2 The following argument is due to V. Lafforgue. 
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The integers Aj^ M+k can again be evaluated as contour integrals of the generating 
series F k (t). In analogy with eq. (|7.13|) , the expressions (|7.24j) , (|7.25|) yield 

f dt 1 j dp M-Np 

At 



h u -''' ~ I 2nit K + 1 p k ~ J 2mp KM + k + l (1 - p) K +' 
I dt 1 r dp M - Np 
A KM+ k = f 2nitK+1 f 27r i p KM +k+ 2 (i ( 7 - 28 > 

(remember that M = N — 1), and finally 

k(KN + k-l)\ + (k + l)(KN + k)\ 

KM+k K\{KM + k)\' KM+k K\{KM + k + l)V 1 } 

The corresponding amplitudes B^ M+k read asymptotically 

B KM+ k « B-{k/M) = k(M + iy k ' M (k = l,..., M), 

(7.30) 

B£ M+k *> B+(k/M) = (k+ (* = 0,...,M-1). 

These expressions agree with eq. (|6.8D , since k/M = Fiac ± (k/p c ) in both cases. We have 
in particular = M/iV = p c , and 5 + (0) = iV/M = l/p c . 

For M large, both amplitude functions B ± (x) exhibit the scaling form 

B ± (x)^xM 1 ~ x (0<x<l, M>1). (7.31) 

This expression takes its maximal value, -B^ ax ~ M/(elnM), for x max = l/(lnM). This 
estimate is only by a factor InM below the upper bound of eq. ( 4.13j ), which reads 
l/q c = N = M + 1. 



7.5 A case study: the slope v = 1/5 or p c = 2/5, i.e., M = 2, 
M = 3 

In this section we study the simplest case of a rational slope which belongs to neither of 
the classes studied before, namely M = 2 and M — 3. This example corresponds to the 
slope v = 1/5. We have N = 5, p c = 2/5, and £ c = 108/3125. We shall apply successively 
the three techniques exposed in sec. ||. 

We start with the approach developed in sec. |5.1| , based on the continuity of the path. 
The functions involved in the linear system (|5.7|) read 




(7.32) 
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The solution of this system leads to the generating series F^(t) of eq. ( |5.4|) in the form 



p ± _ "1 ^22 ~ ^2 ^12 p ± _ ( ~ r 2 Wl ~ "1 ^21 f 7 no\ 

Lt 11 ( - t 22 ~~ ( - J 21 < - T 12 Lr ll ( - ;r 22 — ( - T 21 ( - T 12 



^11^22 ~~ «■ 



The Taylor expansions of these quantities lead to the integers Af, namely 

F+ = 1 + 7t + 99 t 2 + 1768 t 3 + 35530 t A + 766935 1 5 + 17368680 1 6 + 407139120 t 7 

+ 9794689506 t 8 + 240455164510 t 9 + 5999744185435 t 10 + • • • , 
F+ = 2 + 23 1 + 377 t 2 + 7229 1 3 + 151491 1 4 + 3361598 1 5 + 77635093 1 6 



+ 1846620581 t 7 +44930294909 t 8 + 1113015378438 t 9 +27976770344941 1 10 + - 



F- = 1 + 5£ + 66 1 2 + 1144 1 3 + 22610 1 4 + 482885 1 5 + 10855425 1 6 + 253086480 
+ 6063379218 t s + 148365952570 t 9 + 3692150267960 1 10 + ■ ■ ■ , 

F 2 = 2 + 19 1 + 293 1 2 + 5452 1 3 + 1 12227 1 4 + 2460954 1 5 + 56356938 1 6 

+ 1332055265 t 7 + 32251721089 1 8 + 795815587214 1 9 + 19939653287183 1 10 + ■ ■ ■ 

(7.34) 

An arbitrary number of terms can be evaluated with the help of a computer. We have, 
however, found no simple closed form for the coefficients A^, such as eq. ( |7.14j ) for M = 1 
or eq. (|7729| ) for M = 1. 



We now turn to the algebraic trick of sec. |5.2| . For < t < t c , the polynomial 
equation ( |5.13| ), i.e., u 2 (l — u) 3 = t, has two small solutions u\ and U2, such that ui < 
< u\ < 2/5. The relation between t and any of the F^ can be obtained by means of 
an algebraic elimination of U\ and 112 between eqs. (5.13) and ( 5.21| ), a tedious task that 
we prefer to leave to the computer. We thus obtain 

= 



= t 6 (F+) 10 - t\F+) 7 - llt 3 (F+) 5 - t 2 (F+) 4 - 7t(F+) 2 + F+ - 1, 

= t n (F 2 + ) 10 + 10t 10 (F+) 9 + 45t 9 (F+) 8 + t 7 (120£ + 1)(F+) 7 

+ 6t 6 (35t + 1)(F+) 6 + t 5 (252t + 17)(F+) 5 + 30t 4 (7t + 1)(F+) 4 

+ 5t 3 (24t + 7)(F+) 3 + t 2 (45t + 26)(F+) 2 + (10t 2 + lit - 1)F+ + t + 2, 

= t 4 (Fr) 10 - 3t 3 (Ff) 8 + 3t 2 (Ff) 6 + llt 2 (Ff ) 5 - t(Ff ) 4 - 4t(Ff ) 3 + Ff - 1, 

= t 9 (F 2 ) w - 9t s (F 2 -) 9 + 36t 7 (F 2 ~) 8 - 84t 6 (F 2 -) 7 + 126t 5 (F 2 ~) 6 

- 2t\t + 63)(F 2 -) 5 + 3t 3 (3t + 28)(F 2 -) 4 - t 2 {17t + 36)(F 2 ~) 3 

+ t{m + 19)(F 2 -) 2 - (9i + l)F 2 + t + 2. (7.35) 

The algebraic relations between F^ and p are slightly more complicated. We assume 
that p > 2/5. The elimination of t and u\ and U2 between eqs. ( |5.13|) , (|5.21|) , and ( |5.9| ) 
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leads to 



= p 3 (F + f - 10p 3 (F + ) 5 + 3p 2 (15p- 1)(F+) 4 -p(5p- l)(25p- 1)(F + ) 3 

+ 3p(15p - l)(5p - 2)(F+) 2 - 10p(5p - 2) 2 F+ + (5p - 2) 3 , 

= (F-) 6 -(10p-9)(F-) 5 + 3(p-l)(15p-ll)(F-) 4 

- (p - l)(125p 2 - 185p + 63)(F~) 3 + 3(p - l) 2 (15p - ll)(5p - 2)(F") 2 

- (p-l) 2 (10p-9)(5p-2) 2 F- + (p-l) 3 (5p-2) 3 . (7.36) 

The above algebraic relations ( |7.35|) , (|7.36| ) have respective degrees deg t (2, 5) = 10 
in the Fjt and deg p (2, 5) = 6 in F ± , in agreement with eq. (|5.22|) . These relations also 
allow an investigation of the critical behavior of quantities as p — > 2/5. It turns out that 
all the critical amplitudes defined in sec. |6| can be expressed as rational functions of one 
single irrational number: 

C = 10 1/3 , (7.37) 

namely 

125 3125 
(Fl+)c = 7^(2 + 2 C-C 2 ) = 1-929723, (F+) c = — (2 - 10C + 5( 2 ) = 5.889270, 

25 625 
(Fr)c = -(C - 1) = 1-603382, (F-) c =— (3 -0 = 4.893318, 

Q £ TQ 1 or 

$+ = = 5.787037, $+ = (4 - 2C + C 2 ) = 29.020380, 

1 108 '2 11664 v s s ) 

125 3125 

$7 = = 3.472222, $2 = ( 2 + = 20.034890, 

36 648 

$+ = -^(10 - 2C + C 2 ) = 2.870200, $~ = - (5 + C) = 2.384812, 

18 3 

1 1 

C + = — W10 - 2C + C 2 ) = 0.994266, C~ = —^(5 + () = 0.826123. (7.38) 

6v3 5v3 



To close up the study of this example, let us consider the duality approach of sec. 1573 



We have lg = {1,3} and I = {1, 2, 4}. We again focus our attention on quantities with 
the superscript +. Eqs. (|5.40|) , (|5.41| ) yield four equations for the four unknown functions 
r 1 , . . . , r 4 . 

Ft = l + t(F+F+ + F 2 + F 3 + ), 

F + 2 = (Ftf + tFpt 
Ft = FtFt, 

Ft = FtFt- (7.39) 

This non-linear system can be solved in closed form. The trick is to define a variable 
x = Ft /(Ft) ■ The functions Ft and the variable t can then be expressed as rational 
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functions of x: 



x 3 (x 



{x 2 + x - 1) 



5 • 



Ft 



F + 



x 2 + x 



X 

„2 



(x 2 + x - iy 

X 

{x 2 + x- l) 3 



3 — 2 ' 

Ft = <^4^. (7-40) 

In this case, the Riemann surface of t and of the F 1: . . . , F~^ has genus 0. It is a 
branched covering of the t-sphere of degree deg t (2, 5) = 10, with x being the uniformizing 
variable. The physical region < t < t c = 108/3125 corresponds to 1 < x < x c = 
(2 + C 2 )/6 = 1.106931, with dt/dx = at x = x c . We close up this section by observing 
that in general the genus of the Riemann surface where all the functions F~l and t are 
uniform grows very rapidly with M and M. 



7.6 Directed scaling limit: v — ► 1, i.e., p c — ► 0, and p — > 

We now investigate the situation where p and p c are simultaneously small. We refer 
to this the 'directed scaling limit', since the random walks are almost perfectly 

directed toward the rightmost limit (x = t) of phase space [see Figure 1]. This non- 
conventional scaling limit turns out to be characterized by non-trivial scaling laws, in the 
three regimes described in section 2. 
We introduce the scaling variable 

£ = -■ (7.41) 

Pc 



We consider first the convergent regime (£ > 1). The explicit results derived in sec. U~2 



for the family of rational slopes with M = 1 exhibit scaling behavior as p c 1, i.e., 
iV 3> 1. Indeed the small root U\ of eq. (|5.13| ) scales as u\ ~ £iPc, where < £i < 1 is 
related to £ > 1 by 

Eqs. ( 1730 ), ( imp then yield 



£ ie - 4l =£e^. (7.42) 

F^^e^ 1 , (7.43) 

F ± ^(Z-£i)Pc, (7.44) 

n«ei(e-e0Pc- (7-45) 
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More generally, the algebraic formalism of sec. 5^2 can be worked out for an arbitrary 
rational slope p c = M/N in the directed scaling limit, i.e., M finite and N ^> 1. The M 
small roots u a of eq. fl5.13| ), with a — 1, ... , M, now scale as u a ~ £, a Pc, where £ Q is the 
solution such that |£ a | < 1 of 

£ Q e - k = w a "^e^, (7.46) 

with the definition (|5.15|) . As a consequence, solving the linear system (|5.21 ) amounts 
to inverting the discrete Fourier-transform matrix Sk,e = u kl with k,£ — 1, . . . , M. We 
have (S' -1 )^ = u~ H /M, so that finally 

1 M 

*f « rf-*«e-«)-* - £ w-^ - 1 ^. (7.47) 

The scaling form (|7.44j ) for the survival probabilities F ± thus holds true for any 
(rational) slope in the directed scaling limit. The probability II still scales as p 2 c1 but 
with a non-trivial coefficient: 



n«(e-ei)(^£6»)rf- (7.48) 



The critical behavior of the survival probabilities as p — > p c , i.e., ^ — > 1, is character- 
ized by the amplitudes 

$ ± « 2. (7.49) 

We thus obtain the scaling behavior of the amplitudes C ± of eq. ( 2.29| ) in the marginal 
regime (p = p c <C 1 or £ = 1): 

C± « (2 Pc ) 1 /2 ; ( 7 . 50) 

in agreement with eq. ( |6.16 ). 

Some of the above results can be alternatively derived by means of the integers A^. 
The bounds ( 4.10|) imply that these numbers exhibit a simple scaling as p c <C 1, namely 

At - (7.51) 



The behavior of the for large k agrees with eq. ( 4.12 ), with a trivial modulation 

b£ « i. 

Eq. (|3.15| ) then yields 



i^- (£-0(V«))p c , (7.52) 

with 

d>{z) = £ (7.53) 
fc>i 
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This functional relation between (j> and z is equivalent to 



(7.54) 



with the condition that <fi —>■ for z — > 0. The identity (|7.54 ), which is given in refs. 0, |10 
can be checked by means of the contour integrals 

In eq. (|7.52 ) we have z = £e - ^ with £ > 1, so that eq. (|7.42|) leads to the identification 
</> = £i. The result ( |7.44[ ) is thus recovered, without recourse to the algebraic trick. 
Its validity for any slope in the directed scaling limit, either rational or not, is thus 
established. 




0.4 / 0.6 

P c /P 



1.0 



Figure 16: Plot of F/p against p c jp for p = 0.1, 0.2, and 0.3, illustrating the scaling 



law (|7.56 ) pertaining to the directed scaling limit, shown as a dashed curve. 

The result (|7.44j) can be recast as the following scaling relation for the survival prob- 
abilities 

/■ - / /• - \ 

(7.56) 



F ± ( f ±n 
~ 1 — exp 

v V Vc , 



which is valid throughout the convergent regime of the directed scaling limit, i.e., when 
p, p c < p, and F^ are simultaneously small. Figure 16 illustrates this relation. The ratio 
F/p is plotted against p c /p. The data for p = 0.3, 0.2, and 0.1 (full curves, from top 
to bottom), already presented in Figure 10, are found to smoothly converge toward the 
scaling law (|7.56|) , shown as a dashed curve. 
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To close up this section, we consider briefly the large-deviation regime (p < p c or 
£ < 1). When p and p c are simultaneously small with p < p c , the survival probabilities 
F ± (n,v) can be estimated by inserting the expression ( [7.51 ) for the integers into 



eq. ( |6.12| ). We obtain after some algebra an asymptotic expression of the form ( |2.21| ), 
where the entropy function scales as S(v) ~ p c (£ — 1 — ln£), as it should, and where the 
periodic amplitudes a ± (x) and b ± (x) read 

^•1— Frac = F(x) £l— EracrF(a;) 

a±{X) W WM' b±(x) * (2^)^(1-^)' (7 ' 57) 

8 Discussion 

We have performed a detailed analysis of the statistics of persistent events in the case 
of the one-dimensional lattice random walk in the presence of an obstacle moving ballis- 
tically with velocity v. Both space and time are discrete, so that the underlying lattice 
structure yields a highly non-trivial dependence on the velocity v of the obstacle, with 
discontinuities at rational values of v, for most of the quantities investigated in this work. 
This is especially the case for the limit survival probabilities F ± in the convergent regime, 
and for the amplitudes C ± characteristic of the marginal regime, respectively shown in 
Figures 10 and 14-15. We have obtained a deep insight into the problem by means of 
the algebraic methods exposed in sec. 5. 

We want to emphasize that the statistics of persistent events can be investigated in a 
broader perspective, even within the realm of one- dimensional random walks, constrained 
to remain on the left of a moving obstacle. Consider a general random walk consisting 
of independent identically distributed steps s n , and a moving obstacle, whose position 
is given by an arbitrary function X(t). The following situations illustrate the richness 
of possible behavior of this seemingly simple system [see refs. |l], |], 11, [L2j for related 
questions]. 

• The steps e n have a distribution whose first two moments are finite, and the obstacle 
moves ballistically: X(t) = vt. The Sparre Andersen formalism of sec. |2| applies 
in this situation. The survival probability is described by the scaling laws ( |2.21| ), 
( |2.29|) , and (|2.32|) , associated with the three regimes described in the introduction. 



The entropy function S, the constant C of the marginal regime, and the limit 
survival probability F in the convergent regime, are non-universal, in the sense that 
they depend on details of the distribution of the steps. The standard large-deviation 
formalism allows to determine S by means of its Legendre transform. The constant 
C and the limit survival probability F are given by the formal expressions ( p. 28 
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and Q2.33 ). No closed-form expressions for these quantities are known in general. 



Let us mention, however, that the second expression of eq. (|2.35| ) is quite general, 
with Var e = (e 2 ) — (e) 2 replacing 1 — v 2 . 

The steps e n have a symmetric distribution whose second moment is finite, and the 
obstacle moves according to X(t) = gt 1 ^ 2 , in units where (e 2 ) = 1. In this case the 
survival probability of the walker in the presence of the obstacle exhibits a scaling 
behavior of the form 

F(n,g) ~ n -° {9) (n>l), (8.1) 

where the exponent 6(g) is a non-trivial universal, continuously decreasing function 
of g, given in terms of a zero of the parabolic cylinder function JETl [T^J. The 
marginal regime is recovered as 6(0) = 1/2. The result ( |3.1|) crosses over to the 
large-deviation regime as g — > — oo, where the exponent diverges as 6(g) ~ g 2 /8, 
while it crosses over to the convergent regime as g — > +oo, where the exponent 
vanishes as 6(g) « g e^ g2 ^ 2 /(Stt) 1 ^ 2 . The exponent 6(g) is related to the persistence 
exponent 6(g) defined in the introduction, for the sake of consistency with refs. |], 
g §, by 6(-g) = 



The steps e n have a symmetric broad (Levy) distribution, with long tails falling 
off as p(e) ~ |ep M- , with < fi < 2, and the position of the obstacle obeys 
an asymptotic long-time behavior of the form X(t) rs gt 1 ^. Then the survival 
probability is expected to exhibit a scaling behavior of the form ( |S.1|) , again with 
a continuously varying exponent 6(g), that is universal if the e n are measured in 
appropriate units, but whose expression is not known in general. 

A simple example in this category corresponds to steps having a Cauchy distri- 
bution: p(e) = 1/ (n(l +e 2 )), with an obstacle moving ballistically, according to 
X(t) = vt. The Sparre Andersen formalism again applies to this situation. The 
stability of the Cauchy law under convolution implies that the one-time distribution 
function reads 

P(n, v) = P(v) = [ V dS =\ + - arctan v ' ( 8 - 2 ) 

independently of n. Eq. ( [2.1| ) then leads to f(z,v) = (1 — z)~ p ( v \ hence 

, s T(n + P(v)) n~°^ , * «. 

F ^=^kvdi K wm (n>>1) ' (8 - 3) 

where 6(v) is again a universal exponent, given by 

6(v) = 1 - P(v) = - - - arctan v. (8.4) 

2 7T 
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The marginal regime is again recovered as 8(0) = 1/2. The result ( |3.4|) crosses 
over to the convergent regime as v — > +00, where the exponent vanishes as 8(v) ~ 
1/(ttv). Contrary to the previous situation, the exponent does not cross over to 
the large-deviation regime as v — > —00, since it admits a finite limit 6(— 00) = 1. 
The exponent 8(v ) is again related to the persistence exponent 8(v ) defined in the 
introduction by 9(— v) = 8{v), yielding in the present case 8(v ) = P{v). 
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A Combinatorial proof of the Sparre Andersen iden- 
tity (03) 



In this Appendix we provide an elementary and self-contained combinatorial proof of the 
Sparre Andersen identity Q2.1|). 

We keep notations consistent with the body of the paper, dropping for the sake of 
simplicity the ± superscript and the dependence on v. Let {e n } n= i^ t ... be independent 
identically distributed random variables, and let xq — 0, X\ = . . ., x n — £i+- • -+e n , . . . 
be their partial sums. For n = 1, 2, ... we denote by P(n) = Prob{x n > 0} the probability 
that the n-th partial sum is non-negative, and by F(n) = Probjxi > 0, . . . , x n > 0} the 
probability that the first n partial sums are non-negative. We set consistently F(0) = 1. 

We want to prove that the generating functions 

P(*) = E /(*) = E F (n)z n (A.l) 

n>l 71 n>0 

are related by the Sparre Andersen identity Q2.1|), i.e., 

f(z)=exp(p{z)). (A.2) 

For each n > 1, we define the event A n and n events for i — 1, . . . , n by 

A n = {£! + ■ ■ ■ +e n > 0}, 

= {e i >0,e i + e i+1 >0,...,e i + --- + e n >0, (A.3) 
e< H H £n + £i > 0, • . • , £i H h e n + £i H h > 0}. 

The first observation is that A n = 1J™ =1 B* ra . It is clear from the definition that B^ C 
A n . To prove the reverse inclusion, assume A n holds, i.e., x n > 0. Let i e {1, . . . ,n} 
be such that is the minimum of {x , . . . , £ n _i} (in case of degeneracy, we take the 
smallest such i). We claim that B^ is realized. Indeed, by definition of m, x m — > 
for m G {i, . . . , n — 1}, and (as x n > 0), x n + (x m — > for m e {0, . . . , i — 1}: this 
covers the definition of B^ . 

We apply the inclusion-exclusion principle and write 

n 

P(n) = Prob{A n } = Y,(~) k+1 Prob{B^n...nB l „ fe }. (A.4) 

k=l 0<h<...<i h <n+l 

Now comes the second observation. Choose k > 1 and < %\ < ... < if. < n + 1. 
Then 
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n . . . n { £i . > o, . . . , £i . + • • • + > 0} 

n...n{e ik >0,...,e ik + --- + e n >0,..., 

£i k + ■ ■ ■ + £n + £i + ■ ■ ■ + > 0}. (A.5) 

Indeed all the inequalities on the right-hand-side are among the defining inequalities for 
the left-hand-side. Moreover, it is obvious that the missing inequalities on the right-hand- 
side are sums of the written ones. We have thus succeeded in decomposing B^fl. . .flB^ as 
a product of k events. These events are independent because they involve distinct steps of 
the random walk. Another nice feature is that the event {£j > 0, . . . , £j. + • • - + £i- +1 -i > 
0} has probability F(ij + i — ij). Indeed, the steps are identically distributed, so that 

exchanging the steps do not change probabilities, and {e^ > 0, . . . , e^H >0} 

has the same probability as {si > 0, . . . , E\ + ■ ■ • + Ei j+1 -i j > 0}. We can thus rewrite 
eq. ( |A.4| ) for P(n) as 



Pin) = E(") fc+1 E F ^ ~ • • • F (ik ~ ik-i)F{h -ik + n). (A.6) 

fc=l 0<ii<...<i fc <n+l 



The last step of the proof is as follows. Multiply eq. (|A.6|) by z n and sum over n 



1,2, .. . The left-hand-side yields zp'(z), where the accent denotes a differentiation. At 
fixed k, the /c-uple sum in the right-hand-side can be performed by taking j\ = — H > 1, 
32 — h — h > 1) • • jk-i — h — h-i > 1, jk = h — ik + n > 1 as independent summation 
indices. The power of z can be recast as z^ 1+ '" + ^ k . Each of the sums over ji, . . ., jk-i 
brings a factor f(z) — 1, while the last sum over j\ yields zf'(z). We are thus left with 
the equation 

z P '(z) = j:(l-f(z)) k - 1 zf'(z), (A.7) 



k>l 

i.e., 



P'W - % (A. 
which yields the identity (|A.2j ) by integration, since p(0) = 0, /(0) = 1. 
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Figure captions 

Figure 1: Configuration space of the problem. 

Figure 2: A walk remaining on the left of the wall, thus contributing to the survival 
probability. 

Figure 3: A walk crossing the wall, thus not contributing to the survival probability. 

Figure 4: Periodic amplitudes a(x) and b(x), for p = 0.3 and v equal to the golden 
slope. 

Figure 5: Labeling of crossing edges by an integer k = 1, 2, 3, equal to the number of 
steps of the walker to the left. 

Figure 6: Combinatorial approach: truncated Pascal triangle. 

Figure 7: The probability flow is characterized by integers A^, living on the crossing 
edges [cf. Figure 5]. In the example, we have A 1 = 1, A 2 = 2, Af = 5, and so on. 

Figure 8: The canonical walk. 

Figure 9: Roots of the polynomial equation (|5.13|) for M = 4, N = 11, and p — 1/2. 

Figure 10: Plot of the survival probabilities F against v, for several values of p, 
indicated on the curves. The apparent v — > 1 limit of F(v) is F~(l) = p [see eq. ( |7.5|) 1, 
whereas the jump at v = 1 to F + (l) = 1 is not visible. 



Figure 11: Combinatorial proof of the identity ( p. 30 ). 



Figure 12: Plot of the periodic amplitude B(x), for the golden slope. 

Figure 13: Plot of the periodic amplitude (3(8), for the golden slope. 

Figure 14: Logarithmic plot of the critical amplitude C(v) against v. 

Figure 15: Plot of y(v), defined in eq. (|6.14j) , against v. 



Figure 16: Plot of F/p against p c /p for p = 0.1, 0.2, and 0.3, illustrating the scaling 
law ([7.56 ) pertaining to the directed scaling limit, shown as a dashed curve. 
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